Hitting times for the stochastic wave equation with 

fractional-colored noise 
(N' 

Q . Jorge Clarke De la Cerda ^'^ * Ciprian A. Tudor ^' ''' 

^ CPMA, Departamento de Ingenieria Matematica, Universidad de Concepcion, 
c^ I Casilla 160-C, Concepcion, Chile. 

jclarke@udec.cl 

^ Laboratoire Paul Painleve, Universite de Lille 1 
F-59655 Villeneuve d'Ascq, France. 
tudor@math.univ-lillel.fr 



(N 



OO 



^ ■ February 27, 2013 

-)— » ' 



Abstract 

We give sharp regularity results for the solution to the stochastic wave equation 
with linear fractional-colored noise. We apply these results in order to establish upper 
and lower bound for the hitting probabilities of the solution in terms of the Hausdorff 
measure and of the Newtonian capacity. 
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- - -' 1 Introduction 

The recent development of the stochastic calculus with respect to the fractional 
Brownian motion (fBm) naturally led to the study of stochastic partial differential equations 
(SPDEs) driven by this Gaussian process. The motivation comes from the wide area of 
applications of the fBm. We refer, among others, to [10], [11], [13], [16] and [TS]. The 
purpose of our paper is to study the stochastic wave equation driven by fractional-colored 
Gaussian noise. Our work is situated somehow in the continuation of the line of research 
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which concerns SPDEs driven by the fBm but in the same time it follows the 1 research 
initiated by Dalang in [4] which treats equations with white noise in time and correlated in 
space. More precisely, we consider a system of k stochastic wave equations 

-^{t,x) = Aui{t,x) + Wi{t,x), tG [0,r],xGM'^ (1) 

with initial condition Ui{t, x) = and ^(0, x) = for every x € M'^ and for every i = I, ..., k. 
The driving Gaussian process behaves as a fractional Brownian motion in time and has 
spatial covariance given by the Riesz kernel. More precisely 

E{Wi{t,A)Wj{s,B)) = 6ijRH{t,s) [ [ f{x-y)dxdy 

JaJb 

for every t,s £ [0, T] and A, B Borel sets in M where / : M — )• M_|_ is the Fourier transform 
of a non-negative tempered measure /i on W^ whose density with respect to the Lebesque 
measure is |^|~' '^^ < /3 < d. Above (Jj.j denotes the Kronecker symbol. 

The equation ([T]) has been recently studied in [2]. It is has been proven that ([T]) 
admits a unique mild solution if and only if /3 < 2H + 1 which extends the result obtained 
in [3] in the case H = h. The purpose of this work is to analyze further the solution of 
([1]). We will actually give sharp results for the regularity of it, in time and in space, and 
we apply these regularity results to study the hitting probabilities for the solution u to 
([1]). More precisely, given a Borel set ^ C M'^ we want to determine whether the process 
{u{t,x),t € [0, r],x G M.'^) hits the set A with positive probability. Recently, there has been 
several papers on hitting probabilities, and more generally speaking, on potential theory 
for systems of SPDEs. We refer, among others, to [5], [6], [7J, [8J or [I3j. The study 
of hitting probabilities for stochastic partial differential equations with fractional noise in 
time is new. As far as we know, only the paper [15] treated this problem. Actually, in this 
reference the authors give upper and lower bounds for the hitting times of solution to a 
system of stochastic heat equations on the circle with fractional noise in time. 

Our aim is to make a new step in this research direction. As we mentioned before 
we make a potential analysis of the solution to the stochastic wave equation with fractional- 
colored noise. That means, the noise behaves as the fractional Brownian motion with respect 
to the time variable and it is a "colored" non-white spatial covariance. In our work this 
spatial covariance will be described by the Riesz kernel. It is know classical the fact that 
in order to obtain results on the hitting times of a stochastic process, a detailed analysis 
of the behavior of the increments of the process is needed. We address this question in 
our paper and we find the following: the solution u{t,x),t € [0,r],x G M"^ to ([TJ is Holder 
continuous of order 2H + 1 — /3, /3 G {2H — l,d A 2H + 1) in time as well as with respect to 
the space variable. This generalizes the result obtained in [8] and [9] for the wave equation 
with white noise in time and Riesz covariance in space. Although the main lines of our 
work follows the approach of [9], we stress that, as usually, the fractional cases involves 
more complex calculation and the techniques used in the standard white noise case need to 
be substantially adapted, this is mainly due to the nature of the noise and to the structure 



of the Gaussian space associated to the noise. We will point out later in our paper, how 
the fractional noise involves more complexity than in e.g. [6j or [9j. Moreover, the study 
of the solution to the wave equation is generally recognized to be more difficult that e.g. 
the solution to the heat equation, due to the appearance of the trigonometric functions and 
this is also the case in our work. 

We mention that there are more or less general criteria to determine the hitting 
times for a stochastic process. Such criteria have been given in [3], [6], [7] or [9] among 
others. We will use the approach in [3] because it concerns Gaussian processes and fits well 
with our context (note that the solution to ([T]) is Gaussian). 

Our paper is structured as follows. Section 2 contains some preliminaries, we briefly 
describe the basic properties of the Gaussian noise and its associated Hilbert space, we list 
the elements of the potential theory that we will use in our paper and we will recall some 
fact related to the solution to the stochastic wave equation with fractional-colored noise. In 
Section 3 we analyze the Holder regularity of the solution with respect to its time and space 
variables. Section 4 is devoted to the study of the hitting probabilities for this solution, 
based on a criterium in [3]. 

2 Preliminaries 

This section is devoted to introduce the basic notion that we will need throughout the 
paper. We first introduce the canonical Hilbert space associated to the fractional-colored 
Gaussian noise. In the second part we present the basic elements related to the potential 
theory that intervene in the last section. 

2.1 The canonical Hilbert space 

We denote by C^(\R. ~^^) the space of infinitely differentiable functions on R ^-"^ with 
compact support, and 5(]R'^) the Schwartz space of rapidly decreasing C°° functions in M'^. 
For if G L^(M ), we let Tip be the Fourier transform of (p: 



We begin by introducing the framework of [4j. Let ;U be a non-negative tempered 
measure on W^, i.e. a non-negative measure which satisfies: 

/ 1 l—JTv? ) /^(^O < oo> for some / > 0. (2) 

Since the integrand is non- increasing in /, we may assume that / > 1 is an integer. 
Note that 1 + |^|^ behaves as a constant around 0, and as |^p at oo, and hence ([2]) is 
equivalent to: 

/ /u(c?0 < oOi and / /^('^OttttT < ^xd, for some integer I > 1. (3) 

^l€l<i ^l€l>i 1^1 



Let / : M — )• M+ be the Fourier transform of /i in 5'(M ), i.e. 
/ f{xMx)dx = f T^iOKdO, Vy. G cS(R'^). 
Simple properties of the Fourier transform show that for any ip,Tp £ 5(M ) 



ip{x)f{x - y)il){y)dxdy = / T'f{OTip{i)^i{di). 

An approximation argument shows that the previous equality also holds for indicator 
functions (p = 1a, ip = 1_b, with A, B £ ^^^(M'^), where 13b{W^) is the class of bounded Borel 
sets of M.'^: 

f [ f{x-y)dxdy= [ TUiO'^^^i^f^idO- (4) 

Now we introduce the fractional Brownian motion (fBm) with Hurst index H G 
(0, 1). This is a zero-mean Gaussian process (-B/^)te[o,r] with covariance 

R^it, s) := ^ (t^^ + s2^ -\t- sp^) , t,se [0, T]. 

Let us denote by Ti the canonical Hilbert space associated with this Gaussian pro- 
cess. This canonical Hilbert space is defined as the closure of the linear space generated by 
the indicator functions l[o,t]>i G [Oi^] with respect to the inner product 

(l[0,t])l[0,s])w = RH{t,s). 
It is well known that for H > 1/2 we have the expression 

RH{t,s) = an / / |u — up ~'^dudv (5) 

Jo Jo 

for every s,t £ [0,T] with an '■= H{2H — 1). More generally, for H > 1/2 and every 
ip,(t)e'H = n ([0, T]) we have 

{iIj,4>)h = (^H / / i^{u)(f){v)\u — v\ ^ dudv (6) 

JO Jo 

As in [Ij, on a complete probability space (i7, J-", P), we consider a zero-mean Gaus- 
sian process W = {Wt{A);t >0,A£ Bh{R'^)} with covariance: 

B{Wt{A)Ws{B)) = Rnit, s) [ [ f{x - y)dxdy =: {l[o,t]xA, 1[o,s]xb)hv. (7) 

JAJB 

Let £ be the set of linear combinations of elementary functions l[o,t]xyi) t > 0,A £ 
Bb(M'^), and HV be the Hilbert space defined as the closure of £ with respect to the inner 



w 



product (•, ■)n'P- (Alternatively, T-LV can be defined as the completion of Cq°{W^^^), witli 
respect to the inner product (•, •)nVi see [1].) 

The map lro,tixA '~^ ^t(^) is an isometry between £ and the Gaussian space H 
of W, which can be extended to TiV. We denote this extension by: 

ip^W{(p)= / ip{t,x)W{dt,dx). 

Jo JR'i 

In the present work, we assume that H > 1/2. Hence, ([5]) holds. From ([4]) and 
it follows that for any (/?, ■0 € £, 

{¥^,i')HV = OiH I I I I ip{u,x)TJ){v,y)f{x -y)\u-v\^^''^dx dy du dv 
Jo Jo Jw Jr'' 

roo foo c 

\2H-2, 



= an / / T^{u,-){i)F^lj{v,-){i)\u-v\^^-^^i{di) dudv. 

Jo Jo Jm.'^ 

Moreover, we can interchange the order of the integrals dudv and fJ-{d^), since for 
indicator functions if and ip, the integrand is a product of a function of {u, v) and a function 
of ^. Hence, for ip,ijj (^ £, we have: 



{ip,^)HV = aH / / T^{u,-mT^lj{v,-m\u-v\'"-'dudv ii{di). (8) 

JR'* Jo Jo 

The space 1-iV may contain distributions, but contains the space \l-iV\ of measurable 
functions (/? : M+ x M'^ ^^ M such that 

|2 . „ / / / / !,./„. ™MI,.^„. „.MJ^/_ „.M„. „.\2H~2, 



\W\\\HV\ '■= '^H j III W{u,x)\\ip{v,y)\f{x -y)\u-v\ dx dy du dv < oo. 
Jo Jo JR<* JR'i 

2.2 Elements of the potential theory 

Our aim is to analyze the probability 

P {u{i) n A) / 

where u is the solution to ([1]), / is a Borel set included in [0,T] x M'^ and yl is a Borel set 
in M . Here u[I) means the image of / under the random map (t, x) — )• u{t, x). 

We will briefly present the notion of the potential theory that we will need in our 
paper. For all Borel sets F C M we define V{F) to be the set of all probability measures 
with compact support included in F. For all ^ € V{W^), let us denote by //^(/i) the so-called 
/3-energy of the measure /U defined by 



Ip{^l)= / / Kp{\\x-y\\)^l{dx)^l{dy) (9) 



where 



K^{r) = { 



r-^ if /3 > 0; 

log{^) if/3 = 0; 
1 if /? < 0. 



(10) 



Here Nq is a constant. 



For all /3 € M and F € S(M ) we define the /3-diniensional capacity of F by 

-\ -1 



Cap«(F) 



inf I^(/i) 



(11) 



with the convention l/oo := 0. The /3-diniensional Hausdorff measure of the set F ^ B 
is given by 



riRiF) = lim inf 



(12) 



J^(2r,)^; F C U S(x„ri), supr, < e 

where B{x,r) denotes the Euclidean ball of radius r > centered at x € M . When /3 < 
the /^-dimensional Hausdorff measure of F is infinite by definition. 



2.3 The stochastic wave equation with hnear fractional-colored noise 

Consider the linear stochastic wave equation driven by an infinite-dimensional fractional 
Brownian motion W with Hurst parameter H € (0, 1). That is 

— ^(t,x) = Au(t,x) + W(t,x), t>0,xeM.'^ (13) 

u{0,x) = 0, xGM'^ 
-^(0,x) = 0, xeM.'^. 

Here A is the Laplacian on R'^ and W = {Wt{A);t > 0,A e Bb{R'^)} is a centered Gaussian 
field with covariance 

B{Wt{A)Ws{B)) = RH{t,s) [ f fix- y)dxdy, 

where Rh is the covariance of the fractional Brownian motion and / is the Riesz kernel. 

Let Gi be the fundamental solution of uu — Aw = 0. It is known that Gi{t, •) is a 
distribution in 5'(M ) with rapid decrease, and 



j'Giit,-m 



sin(tlel) 
1^1 ' 



(14) 



for any ^ G M*^,* > 0, d > 1 (see e.g. [19j). In particular, 

Gi{t,x) = -l{|x|<i}, ifd=l 

Giit,x) = , l{|x-|<t|, ifd = 2 



27r^ 



\x 



Gi{t,x) = Cd-dt, lid = 3, 

where at denotes the surface measure on the 3-dimensional sphere of radius t. 

The solution of (fT3]) is a square-integrable process u = {u{t,x);t > 0, x E M^} 
defined by: 

u{t,x)= / Gi{t-s,x-y)W{ds,dy). (15) 

Jo JR'i 

By definition, u{t, x) exists if and only if the stochastic integral above is well-defined, 
i.e. gtx '■= Gi{t — •,x — •) G l-LV. In this case, E|M(i,x)p = ||5ix|||^p- 
The following result has been proved in [2]. 

Theorem 1 The stochastic wave equation [T3\) admits an unique mild solution {u(t, a;))(gro_r]^. 
if and only if 



1 \^+3 



,1 + iei 

Remark 1 Note that Iil6\) is equivalent to: 



^i{dC) < oo. (16) 



/ fj.{d^) < oo, and / ^dS.) T-T^rrrT < oo. (17) 

As mentioned in the Introduction, we will consider throughout the paper that the 
spatial covariance of the noise W is given by the Riezs kernel. That means the measure ii 
is 

d^{^) = ICr'^+^dC with /3 e (0,d). 

In this case the kernel / is given by 

/(^) = |^|-/3with/3G(0,d). 

Note that in the case of the Riesz kernel, condition (|16p is equivalent to 

/3g (0,dA(2F + l)). (18) 

Remark 2 Since H > -^ and so 2H + 1 G (2, 3), for dimension d = 1, 2 we have (3 G (0, d) 
while for d > 3 we have /? G (0, 2H + 1). 



3 Regularity of the solution 

3.1 Time regularity 

In this part we will focus our attention on the behavior of the increments of the 
solution u{t, x) with respect to the variable t. We will give upper and lower bounds for the 
L^-norm of this increment. Usually, obtaining upper bounds is recognized to be easier than 
obtaining lower bounds, this is also the case in our work. Actually, to get the sharpness 
of the regularity of u with respect to the time variable, we need to impose a stronger 
assumption than psp on the parameters /3 and H (condition ()19p below). This is due to 
the characteristics of the scalar product in the TiV. 

We will start with the following useful lemma that gives an explicit expression for 
the Ti norm of the cosinus and sinus functions. These norms will widely appear further in 
our computations. 

Lemma 1 Let f{x) = cos(x) and g{x) = sinx for x G M. Then for every a, 6 G K, a < b 

fb—a j-b—a 

ll/l(a,f<)ll|{ = CK-ff / dvcos{v)v'^ ~^(6 — a — f) + a// cos(a + 6) / dvv^ "^ sin(6 — a — u). 

JQ Jo 



and 

pb—a j-b—a 

\\g'^(a,b)\\H~ '^H / dvcos{v)v ~ {b — a — v) — an cos{a + b) / 
Jo Jo 



dvv sin(5 — a — f) 



Proof: Using the expression of the scalar product in the Hilbert space % and the trigono- 
metric identities 

cos(n it f ) = cos n cos f =F sin u sin w and 
sin{x) — sin{y) = 2cos{ — - — )sin{ — - — ) 

we can write 

pb pb 

\\f'^(a,b)\\n ~^ \\s'^(a,b)\\'H = (^H / W — v\ ~ {cOS U COS V -\- siu U siu v) dudv 

J a J a 
/•b rb 

= an du dv\u — v]"^ ~'^ cos{u — v) 

J a J a 

pb pu—a 

= 2aH du dvcos{v)v'^ ~^ 

Ja Jo 

pb—a 

= 2aH / dvcos{v)v ^ {b — a — v) 
Jo 



where we made the change of variables v = u — v in the integral dv above and we computed 
the integral du. Similarly 

/•b j-b 

\\f^{a,b)\\'H.~ \\9^(a,b)\\'H. ~ ^H du dv\u — v\ ~ (cos ti COS f — sin u sin f ) 

J a J a 



fb j-b 

OLH I W ~ v\ ^ cos(ti + v)dudv 



a J a 



and by the change of variable v = u — v m. the integral dv, 

pb pu—a 

\\f'^{a,b)\\H-\\9'^{a,b)\\H = "^oiH du dv cos{2u - v)v^^~'^ 

Ja Jo 

rb—a fb 



po—a po 

2aH / dvv ^ / ducos{2u — v) 

Jo Jv+a 

pb—a 

an / dw^^-"^ (sin(26 - v) - sin(2a + v)) 
Jo 

pb—a 

2aHCOs{a + b) / dvv ~ sin(6 — a — f). 

^0 



Remark 3 As a consequence of the Lem.ma{l\we deduce the following 

i. For any x > the quantity j^ v^^~'^ cos{v){x — v)dv is positive (it is the sum of two 
norms). 

ii. For every a, 6 G M, a < 6 

pb—a 

ll/l(a,6)llw<2a// / dvcos{vy^~^{h-a-v) 
Jo 

Hi. For every a, 5 £ M, a < b 

pb—a 

ll/l{a,fe)l|-H — 2a_ff cos(a + 6) / dvv ^ sm{b — a — v). 

Jo 

Later, we use also the following lemma. 

Lemma 2 For every a, 6 G M with a <b, 

b fb 

dudv s\'n.{u — v)\u — v\ ~ = 0. 



Proof: This follows from the trivial equality 

b fb fb rb 

sin[u) cos{v)\u — v\ ^ dudv = I I sni{v) cos{u)\u — v\ ^ dudv. 



Concretely, we will prove the following result concerning the time regularity of the 
solution to p^ . We mention that, in the rest of our paper, c,C... will denote generic 
positive constants that may change from line to line. 

Proposition 1 Assume that 

P ei2H -l,dA{2H + l)). (19) 

Let to,M > and fix x G [—M,M]'^. Then there exists a positive constants ci,C2 such that 
for every s,t £ [to,T] 



ci 



|t _ 5|2^+l-/3 < E \u{t, X) - U{S, X)P <C2\t- S|2^+^-^. 



Proof: Let h > and let us estimate the L'^{Q)-novui of the increment u{t + h, x) — u{t, x). 
Splitting the interval [0, t + h] into the intervals [0, t] and [t, t + h], and using the inequality 
\a + 6p < 2(a2 + b'^), we obtain: 

B\u{t + h,x) -u{t,x)\^ < 2{\\{gt+h,x - 9t,x)l[o,t]\\Hr + \\9t+h,x^t,t+h]\\Hr} 

=: 2[E,4h) + E2{h)]. (20) 

The first summand can be handled in the following way. 

Ei^t{h) = an fJ'idO / dvdv\u - v\'^^''^J='{gt+h,x - gtx){u,-){0 

jRd Jo Jo 



xl'{gt+h,x - 9tx)iv, ■){0 



an I KdO I I ' 

Jw^ Jo Jo 



dudv\u - vY"-^[TGi{u + h, 0(0 - ^Gi{u, ■){C}\ 



xTGi{v + h,-m-TGi{v,-m 

ft ft 
an / / dudv\u — v\ Ih, 



Jo Jo 
where 



h = [ fiidO[TGiiu+h,-m-j'Giiu,-m][j'Giiv+h,-m)-TGi{v,-m] 

jR'i 

(sin((n + km) - sin(nlel)) (sin((7; + km) - sm{vm 



KdO- 



10 



Using the last trigonometric identity presented before we obtain 



t ft 



\2H-2 



KdO 



M'^? ,(2u+hm, A2v+hm- 



lep 



COS 



de 



Ei^t{h) = an / / dudv\u — v 
Jo Jo 

and by making the change of variables u = {2u + h)\^\,v = {2v + /i)|^| 



-) cos(- 



sm{h\S,\f cos((2u + h)\S,\) cos((2i; + h)\^\), 



El Ah) 



c ■ an 



I 



sin(/i|^|) / / dudv\u — v\ cosucost; 



\^\d-li+2H+2 



m 



h\^\ 



\^\d-l3+2H+2 

and using Lemma [H 
Ei,t{h) = 



sin(/i|^|) ||cos(-)l(ft|g|,(2i+ft)|?|)(-)ll-H' 



(21) 



d^ 



\^\d-l3+2H+2 



sin(/i|e|)2 X 



r.2t|e| 



2tiei 



cos(u)t;^^"^(2i|^| -v)dv 



+ cos(2i|e| + 2h\^\) / t;2/f-2(g.^(2t|C| - v)) 
Jo 



d^ 



■sm{h\^\f X 



- sin(2i|e|)(2t|e|)2-f^-i + (2/? - 1) / sm{v)v'^^-^dv 

Jo 

+ cos(2t|e| + 2/i|e|) / v^"-\sm{2t\^\ - v)) 
Jo 

where we use integration by parts. By Remark [31 point ii. we have the upper bound 

dC 



(22) 



El,t{h) < C-OH 

Jw 



l^ld-l3+2H+2 



sm{h\(\y 



2*51 



2i|^| / cos(f)f^^-^df-sin(2t|^[)(2t|^|) 





2H-1 



r2t\l;\ 

+{2H-l) / siu{v)v'^^-'^dv 
Jo 



11 



We will treat separately the three summands above. Concerning the first one, 






\^\d-l3+2H+2 





'^^ sin(|C|)2 / ' cosivy-^dv 



\^\d-l3+2H+l 
dC 










2t|S| 

h 



cos(u)w dv 



using condition p8|) and the fact that the integral J^ cos{v)v'^^^'^dv is convergent (this 
implies that the function x € [0, cxd) -^ J^ cos{v)v "^dv admits a limit at infinity and it is 
therefore bounded). On the other hand 



/ 



di 



cth 



Uld-l3+2H+2 
3-/3 f d^ 



smih\C\f sm{2t\^\){2t\^\f^-^ = ah^-^ 



l?l<i 1^1 



d-/3+3 



sindeD^sinf^l+Q/i^-/' 



sm(l^l) sml-^ 

2tiei 



i5i>i 



|^|d-/3+3 

d^ 

|^|a!-/3+3 



sin(l^l)^ sin 



/i 



The second part over the region |^| > 1 is bounded by ch^~" simply by majorizing sinus by 
one. The second integral has a singularity for |^| close to zero. Using that sin{x) < x for 
all x > 0, we will bound it above by 



""L^""''^"'""!^ 



< Cth 



3-/3 



dC 



ICI<i 1^1 



d-/3+3 



^-'^'''-'-'Lw 



d^ 



sm 



2fiei 



2-2H 



sm 



2m 



2H-1 



where we bounded 



sm 



2t|?l 



51: 

2-2H 



^/3+2H-l 



by ctim 



'l\2-2H 



and 



sm 



2m 



2H-1 



by 1. The last 



12 



integral is finite since j3 > 2H — 1 (assumption (|19|) ). Finally 

sin(t>)f ~ dv 



I 



d^ Z"^*!' 



h 



2H+2-I3 



di 



\^\d-l3+2H+2 

,2H+2-l3 I d^ 

|«|<1 lel'^-^+2^+2 



2t|«| 



,.2H-2, 



sin(l^l)^ / sm{v)v^''-^dv 





< h^H+2^^ 
^h2H+2- 



di 



€I>1 



|^|d-^+2fl'+2 



de 



2t|g 

sin(le|)2 / ' sin(t.)t;2^-2(ii; 
sin(lel)^ / " sin(7;)z;2^-2^w 

2t|CI 



|5|<1 |e|'^-/3+2/^+2 
d-^+2_ff+2 



lep 



I sin wit' dw 



l^l>i 1^1 



2t|{ 



sin(w)u du 



(23) 



Using again the fact that J^ sm{v)v'^^~'^dv is convergent is easy to see that the integral 
over the region [^| > 1 is bounded by cj/i^^^^^^. For the intgral over |^| < 1 we make the 
change of variables v = ^ and we get 



.3-/3 



di 



,3-/3 



di 



2t 



sm 



v\i\ 



5I< 



1 lei^-^+i JO 



2t 



sm 



v\i\ 



v^^'-^dv 

2-2H 



sm 



v\i\ 



2H-1 



v^^'-'dv 



< cth 



2/1+1- 



di 



Cl<i 



|^|d-/3+2H-l 



where we have made the same considerations as for the second summand in the decompo- 
sition of Ei^tih). In this way, we obtained the upper bound for the summand Ei^t{h) in 

(EOD 

Ei,t{h) < Ch^"+^-^. (24) 

We study now the term E2{h) in ([20]) (its notation E2{h) instead of i?2,t(^) is due to 
the fact that it does not depend on t, see below). Using successively the change of variables 
u = ^,v = f in the integral dudv and ^ = /i^ in the integral d^, the summand E2{h) can 
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be written as 

rt+h ft+h 



E2{h) = an I I / FGi{t + h - u,-){^)FGi{t + h - v,-){S)\u- v\^^''^du dv ii{dO 

JR<* Jt Jt 



KdO 



fh rh 



lei 



pn pn 

/ / sm(n[^|) sin(f |^|)|u — u 
Jo Jo 



2H-2, 



OiH / I , .„ / / sinfnt^l) sinff |^|)|u — ul dudv 



\i? JO JO 



anh'"'^''-^ I _4§^ / / sm{u\^\)sm{vmu-v\'''-''dudv. 



anh^^ I nrrl^ / / sm{u\^\h) sin {v\C\h)\u-v\'^"-^ dudv 
i{dt 

Let us use the following notation: 

^*(0 = ^ f^ f^^Hu\i\)^Hv\mn-vr-''dudv, t G [0,r],C eM-^. (25) 

By Proposition 3.7 in |2j the term 

iVi(^) = ^ / / sin{u\i\)s\n{v\i\)\u-v\^"-^dudv 

\i\ Jo Jo 

satisfies the inequality 

^ N H+l/2 



NiiO < Ch 



1 + ICI 



with Ch a positive constant not depending on h. Consequently the term E2{h) is bounded 
by 

and this is clearly finite due to (jlSp . Relations (j24p and (|26p give the first part of the 
conclusion. 

Let us analyze now the lower bound of the increments of u{t, x) with respect to the 
variable t. Let /i > 0,j; € [-M^MY^ and t G [to,T] such that t + h£ [to,T]. From the 
decomposition 

E\u{t + h,x) -u{t,x)\^ = II {gt+h,x - 9t,x) l[o,t]\\nr + \\9t+h,x'^[t,t+h]\\nr 

+'^{{9t+h,x - 9t,x) ^0,1], 9t+h,x'^[t,t+h])HV 
we immediately obtain, since the second summand in the right-hand side is positive, 

E\u{t + h,x) -n(t,x)|^ > II {gt+h,x - 9t,x) l[o,t]\\nv + '^{{9t+h,x - 9t,x)\o,t\,9t+h,x\t,t+h])HV 

:= Ei,t{h) + Es,t{h). 
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We can assume, without any loss of the generahty, that t = ^^ Denote E^ i{h) := Ei{h). 
We first prove that 

Ei{h) > c/i^^+i-^ - c'/i2^+2-/3_ (27) 

for h smah enough. Recall that we have an exact expression for Ei{h) (see (j22p ). Actually, 



Eiih) 



d^ 



an 



an 



U\d-(5+2H+2 



Uld-l3+2H+2 
di 



sin(/i|^|)2||cos(-)l(ft|5|,/i|^|+|g|)|||^ 

sin(/i|^|) / / dudv\u — v\ cosucosu 

Jh\^\ Jh\e,\ 

r\S.\ /-I?! 
s\n{h\i\f / / dudvcos{u + h\i\)cos{v + h\i\)\u-v\^"-'^. 
Jo Jo 



|^|a!-/3+2//+2 

By the trigonomtric formula cos(x + y) = cos(2;) cos(y) — sin(x) sin(y) we can write 



Ei{h) 



lei 



d-l5+2H+2 



Mhm 



i?i 



cos(/i|^|)^ / / dudv COS u COS v\u — v\'^ ^ 
Jo Jo 



2H-2 



— 2sin(/i|,^|) cos(/i[^|) / / dudv sin u cos v\u — V 
Jo Jo 

+ sin(/i|^|) / / dudv sin u sin v\u — v\'^ ~^ 
^0 Jo 

:= A + B + C. 

We will neglect the first term since it is positive. We will bound the second one 
above by c/i^^^^^^ . We have (we use Lemma [2] at the third line below) 



sin(/i|,^|) cos(/i|.^|) / / dudv sin u cos v\u — V 
Jo Jo 

H5I /-ICI 



2H-2 



< |sin(2/i|^[)| 



|sin(2/i[C|)| 



c|sin(2/i|^|)| 



dudv sin u cos v\u — v 



2H-2 



Jo 



dudv sin(u + v)\u — v 



\2H-2 



JO 



2H~2 



v''''~'\cos{v) - cos(2|^| - v))dv 
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and thus 
B < c 

< c 



d^ 



|^|d-/3+2H+2 

dC 

Uld-l3+2H+2 



m 



sm(/i|^|)2|sin(2/i|^|)| 

sin(/i|C|)'|sin(2/i|C|)sm(!e|)| 
sm(/i|^|)2|sin(2/i|^I)sm(!e|)| 



2H-2 



z;^''"^(cos(u) - cos(2!^| - v))dv 



/ sm{v + \^\)v'^^''^dv 
Jo 



sin([^|) / cosvf "df + cos(|^[) / sinvv dv\ 



I5l<i 



dj 

\^\d-li+2H+2 



sin(/i|e|)^|sin(2/i|C|)sin(|^|)| 



/ r\<i\ r\i\ \ 

sm(l^l) / cosvv^^ '^dv + cos{\i\) I smvv^^~'^dv\ 



+ c 



/ 



d^ 

\^\d-(i+2H+2 



sm(/.|e|)2|sm(2/.|e|)sin(|^i)i 



sm([^|) / cost;f^ '^ dv + cos{\^\) I svuvv^ ^dv\ 
The part over the set |^| < 1 is bounded by ch^ by shnply majorizing sin(x) by x and 



sm 



(iei)/c 



1^1 ^^o,„,2H-2 



cos VV 



dv + cos{\^\) j^ sm.vv^^ "^dv 



by a constant. Concerning the part 
over the region I .^ I > 1 we bound again ( sin(|,^|) J^^' cosfw^^^^dt' + cos(|^|) Jq siniJu^^^^dt; 
by a constant and we use the change of variables ^ = h^. This part wih by bounded by 

2H+2~P f dS, 



h' 



\>h 



< h^H+2-P 



Uld-l3+2H+2 

dj 

\^\d-P+2H+2 



sin(|e|)^|sin(2|e|)| 



sin(|e|)2|sin(2|e|)| 



since the last integral is convergent at infinity by bounded sinus by one and at zero by 
bounding sin(x) by x and using the assumption /3 > 2H — 1. Therefore 



B < c/l2^+2-/3_ 



(28) 



We bound now the summand C below. In this summand the H norm of the sinus 
function appear and this has been analyzed in [2]. We have, after the change of variables 
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U = TFT,V 



M'^-M' 






C = I _T:^T^^sin{h\^\)^ I I sinmsm{vmu-vf^-^dudv 



di '■I '■1 



> / ^^^'^^(^1^1)'/ I ^Hm^Hy\mu-v?''-^dudv. 



We will use Proposition 3.8 in \^ (more precisely, we will use the inequality two lines before 
the formula (34) in that paper with A; = 0; we notice that the term sin(/i|^|)^ does not 
appear in this proof but by analyzing the step of the proof we can see that this term can 
be added without problems). We will have that, for h small, 

C > I -^^sm{h\i\f-^ I I sm{\i\)sm{v\i\)\u-v\^^-'^dudv 



|^p3^sin(/i|e|)^^2lm 



> / T7±sMhm 



h 



2H+1-I3 f d^ .^r.(h\>:\\^ 



r2HTT^™(^l^l)^ 



= c/i^^+i-^. (29) 



Relations dSH]) and ([29]) imply ^. Now, from relation ^, for every to < s < t < T 
with s,t close enough 

E,it-S) > c(t-s)2^+l-^-c'(t-s)2^+2-/3>£(^_^^2H+l-/3 

if |t — s| < ^. To extend the above inequality to arbitrary values of |t — s|, we proceed as in 
P, proof of Proposition 4.1. Notice that the function g((t, s, x, y) :='E\u{t,x) — u{s,x)\ is 
positive and continuous with respect to all its arguments and therefore it is bounded below 
on the set {{t,s,x,y) G [to,T]'^ x [— M, M]^"^; \t — s\ > e} by a constant depending on e > 0. 
Hence for |t — s| > ^ it also holds that 

On the other side, from (j2ip and (I26p and Cauchy-Schwarz inequality, we obtain 
E3,tW = {{9t+h,x - gt,x) '^[o,t], 9t+h,x'i^[t,t+h])nv 

< II {9t+h,x - gt,x) \o,t]\\Hv\\9t+h,x\t,t+h]\\'HV 

2H+1-I3 . 2H+2-I3 

< Ch 2 + 2 
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Consequently, 

E|u(i + /i,x)-n(t,x)|2 > Ch^"+^-^ -C'h-^^+^^^ 

and this implies that for every s, t G [to, T] and x G [— M, M]'^ 

1 
EKt,x)-n(5,x)|2 > ^|t-s|2^+i-^ if |t-,|<(^-^y. 

Similarly as above, the previous inequality can be extended to arbitrary values of s,i G 

■ 
Proposition ([T|) implies the following Holder property for the solution to (|13p . 
Corollary 1 Assume hl9^) . Then for every x G M the application 

t —?■ u{t,x) 
is almost surely Holder continuous of order 5 G ( 0, 2~^ ) ■ 



Proof: This is consequence of the relations (I2ip and (j26l) in the proof of Proposition [T] 
and of the fact that u is Gaussian. ■ 

Let us make some comments on the result in Proposition [TJ 

Remark 4 

• Following the proof of Theorem 5.1 in 18] we can show that the mapping t — )• u{t,x) is 

not Holder continuous of order — \^ ~^ . 

• When H = ^ the solution of the wave equation with fractional noise in time has the 

same regularity in time as the solution of the wave equation with white noise in time 

(see m, m)- 

3.2 Space regularity 

Let us discuss the behavior of the solution u to the equation (fT3]) with respect to 
the spatial variable. We have 

Proposition 2 Assume il9\) . fix M > and t G [to, T]. Then there exist positive constants 
C3, C4 such that for any x,y £ [—M, M^ 

C3\x - y\^^+^-^ < E \u{t, x) - u{t, y)f < c^lx - yp^+^"'^. 
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Proof: Let z G M*^. We compute 

E\u{t,x + z) -u(t,x)\'^ = \\gt,x+z - gt,x\\np 



t rt 



an I / J='igt,x+z- gt,x)iu,-)iOJ^i9t,x+z - gt,x)iv,-){S,)\u-v\ du dv fi{di) 



= aH [ f \u-v\^^-^dudv f \e-'^<^+''> -e-'^-^\^TGi{u,-)iOJ^Gi{v,-){C) fiidO 
Jo Jo Jr'^ 

= an f f \u - vr-^du dv I |e-«- - i^^^^^ . '^^^{d^ 

Jo Jo JMd l?l l?l 

= :Ei^^{z)+E2,x{z), 

where Ei^x{z) and E2^x{z) are the integrals over the regions \^\ < 1 and \^\ > 1 respectively. 
For the first expresion is easy to see that, using the inequality |1 — e~*^^p < |i^p|zp, we get 
the bound 

Ei,x{z) < C\z\^ [ KdO- 

Jm<i 

Developing the second expresion we get 

Jo Jo J\£,\>1 l?l l?l 

where z ■ ^ means the scalar product in M . Again from Proposition 3.7 in [2] we have that 

/ , Ni/+l/2 

mi) < cuH [j^) 

for any t > 0, |^| > 1, where Nt{S,) is given by (f25]l . Hence, denoting by e 



where we used the change of variables w = ^\z\. This proves the upper bound. 

Let us prove the sharpness of this bound (i.e. the lower bound). We can assume, 
without losing the generality, that f = 1. We note that 

E|u(l,x + z) -u{l,x)\^ >F2{z) 



2aH / / dudv\u-v\'^^ ^ / — 7^3^(1 - cos(^ • z)) 
Jo Jo J\i\>i l?l 



2H-^2 f d^ /I fy- NNsin(ul^l) sin(i;|^|) 
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Condition p7|) implies that 



^,{dO ^ f /i(de) ^ ^_ 



I€I>1 lep -4>llep''"'' 

We apply Proposition 3.8 in [2j (more precisely, the inequality two lines before (34) 
in [2] with k = 0) and we get (note that the result in [2] is stated without the factor 
(1 — cos(.^ • z)) but by analyzing the steps of the proof we can see that this factor may be 
added without problems) 

d^ 



F2{z) > C 1^^^^^ ^7iz^T2lm (1 - ^°«(^ • ^))- 



and by the change of variables ^|z| = w in the integral d^ 

F.(.) > C.— ^/ ^^^^4^(l-cos(u..e). 

J \w\>\z\ \"^\ 

As in the proof of Theorem 5.1 in [8j, we obtain that the integral J|^|>u| i id-0+2H+i (1 ~ 
cos{w ■ e)) is bounded below by a constant. (Notice that /3 > 2H — 1, implies that the first 
integral above is convergent when z is zero, because 1 — cos(x) ~ x^ around zero). Thus, it 
is immediate that 

E|n(l,x + z) - u{l,x)f > Cz2^+i-^ - CV^+2-/3 

and this implies 

E |n(l, x + z)- u{l,x)f > ic7z2H+i-/3 

for \z\ < ^fjT. It is a routine argument to extend the above inequality to arbitrary values of 
\z\ (see e.g. [9], page 22, see also the proof of Proposition [1] before) . 

■ 

We have the following result concerning the Holder continuity in space. We mention 
that it is a little bit more than an extension of Proposition [2j 

Proposition 3 Assume I119\) . Then for any t € [toi^] the application 

X — )• u(t,x) 
is almost surely Holder continuous of order 6 £ I 0, ( — ^ ~ ) '^ ^ ) • 
Proof: We claim that 

E \u{t, X) - U{t, y)f < C\X - y\i^H+l-l3)A2 (3Q) 
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for \x — y\ small enough. From Proposition [51 (j3Up is true when /3 > 2H — 1. When 
/3 S (0, 2H — 1) then it suffices to regards the part of the quantity E \u{t, x + z) — u{t, x)\ 
over the region |^| < 1 (the part over the region |^| > 1 is, as in the proof of Proposition [21 
bounded by cz'^^~^^~^ so by cz^ for z small). It is immediate to see that, using the inequality 
|1 — e~*^^p < iCPl-zp the considered part is less than C[zp f,^,^-.fi{d^). This concludes the 

proof of dSOj)- 

The conclusion is a consequence of Proposition [21 the Gaussianity of u and the 
Kolmogorov continuity theorem. 



Remark 5 

• When H = -^ the solution of the wave equation with fractional noise in time has the 

same regularity in space as the solution of the wave equation with white noise in time 
(see m, m- 

• We distinguish in Proposition 13\ two cases: if f3 ^ (0,2/f — 1) then the solution to U3\) 

has spatial Holder continuity of order 1 (so, it is Lipschitz continuous in the space 
variable) while if P G {2H-l,dA{2H+l)) the Holder exponent is 5 £ (0, ^^+^~^ ) < 1. 

• There is another way to see why the cases j3 G (0, 2H — 1] and (3 £ {2H — 1, d A {2H + 1)) 

need to he separated. Denote by 

gt{z) := E\u{t,x + z) — u(t,x)\ 

and let us study the behavior of gt around z = 0. Let us also assume that d = 1. 
Notice first that gt{0) = and 

n. f* f\ i2ff 2 , , /■ di , ^,sin(u|^|) sin(t;|^|) 

and thus g't^O) = provided that (3 < 2H . Moreover 

II, NO /"*/"*! m-'iA A f ^g I ^, sm{u\C\) sin(^|g|) 

and 



* '■* 2H-2^ J f d^ sin(u|^|) sin(7;|^|) 



ff"(0) = 2aH I I \u-v\^^-'^dudv I 
Jo Jo Jm 

< Ct2aH 



JO JR^ \^\'-^-' lel 1^1 

d^ ( 1 N " ' 5 



which is a finite constant for j3 < 2H — 1. Therefore gt{z) behaves as Cz'^ for z close 
to zero. 
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3.3 Joint regularity 

Let us denote by A the following metric on [0, T] x M*^ 

A {{t, x); {s, y)) = \t- s^^'^^ + \x - yp^+i"^. (31) 

From Propositions [1] and m we obtain the following result: 

Theorem 2 Fix M > and assume /il9\) . For every t,s € [to,T] and x,y £ [—M,M]'^ 
there exist positive constants Ci , C2 such that 

CiA((t,2;);(s,y))<E|u(t,x)-u(s,y)|2<C72A((t,x);(s,y)). 

Proof: The upper bound can be easily obtained by using the upper bound in Propositions 
[U and [2] since 

E |n(t,x) — u(s,y)| < 2'Ei\u{t,x) — u{s,x)\ +2'Ei\u{s,x) — u{s,y)\ 

Concerning the lower bound, it suffices follow the lines of the proof of Lemma 2.1 
in |15) (see also Steps 3 and 4 in the proof of Proposition 4.1 in [9]). We will briefly 
explain the main lines of the proof. The demostration needs to be divided upon three cases: 

\t-s\2H+l-P < c^W_y\2H+l-P U_,\2H+l-P > 4c4 |^ _ |2H+l-/3 ^^^ 4c4 |^ _ |2i/+l-/3 > 
I I — 4c2 \ a \ 5 1 I — ci I ^ I Ci I ^ I — 

\t — s\'^^^^~l^ > j^lx — yl"^^^^^^ with the constants ci, 02,03, C4 appearing in the statements 
of Propositions [1] and [21 The first case can be handled as follows 

E\u{t,x) -u(s,y)\^ > -'E\u{t,x) -u{t,y)\^ -E\u{t,y) -u(s,y)\'^ 

> lc3|x-y|2^+i-/5-C2|t-sp^+i-^ 

> lcs\x-yr+'-^-\c,\x-yr+'-f^ 
= lcs\x-yr+'-f^ 

- 8 ' ^' + 8 C3 ' ' 

> CiA((t,x);(s,y)). 

The other cases follows similarly from Lemma 3.1 in pK], by replacing their expo- 
nents with our exponents. 



Remark 6 The result of Theorem \^ can he stated also in the following form: Fix M > 
and assume \1S^) . For every t,s £ [to,T] and x,y G [—M,MY with (t,x) close enough to 
{s,y), there exist positive constants Ci,C2 such that 

Ci {\t -s\ + \x- y\f''+^~^ < E \u{t,x) - uis,y)\^ < C2 {\t -s\ + \x- y\f''+^-^ . 
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4 Hitting times 

Let us discuss the upper and lower bounds for the hitting probabilities of the solution 
u to equation (J13p . These bounds will be given in terms of the Newtonian capacity and 
the Hausdorff measure of the hit set (see Section 2 for the definition). Let us recall the 
notation: if V = {V{x),x € W^) is a M'^ valued stochastic process then V{S) denote the 
range of the Borel set S under the random mapping x — ?> V{x). 

Our result is based on the following criteria for the hitting probabilities proven in 
[3], Theorem 2.1. 

Theorem 3 Let X = X{t),t G M^ be a M.^ -valued centered Gaussian process and fix 
I C M^. Assume that there exist positive constants 01,02,03,04 such that 



I. 



For every t e I,E [X{t)'^] > oi > 0. 



ii. There exists ai, ..., ajy € (0, 1) such that for every t = {ti, ...jtjy), s = {si, ..., sn) G I 
it holds that 

N N 

02 Y^ \tj - Sjl'^''^ < E \X{t) - X{s)f < 03 ^ \tj - Sjl^"^. 

Hi. For every t = (ti, ...jtjy), s = (si, ..., sjy) G I 

N 

Var{X{t)\X{s)) > 04^] l*i " *jf "'• 
i=i 

Then there exist positive constants 05, og such that for every Borel set A in M'^ 
a^Capk^giA) < P{X{I) n .4 / 0) < aenk^giA) 

where Q = J2f=i ^■ 

Next, we will show that the solution to (|13p satisfies the assumptions of the previous 
result. This will be done via several lemmas. 

Lemma 3 Assume [TPj] and let u he the solution to I73j) . Then for every t € [to,T] and 

xeM.'^ 

Bu{t,xf >C. 
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Proof: Let a^^ be the variance ofu{t,x). We need to give a lower bound for this variance. 
Assume for simplicity t = 1. Then 

(^1 X = E|u(l,x)[ 

= an dudv\u - v\^^'^ j ^ i^.^.^+a sin(^|^|)sm(r;|g|) 

> an dudv\u - v\'^^''^ , „ sin(u|g|) sin(^|^|) 

Jq Jo J\^\<i l?l '^^ 

1 /.I 



pi pi 

Oil 1 Q f/ 

> an sin 1 / / dudv\u — v\ ^ uv = C > 
Jo Jo 



where we used the bound sinx > xsinl for every x G [0,1]. The general case t G \to-,T] 
follows in the same way by doing the change of variables u = j,v = j and then working on 
the domain D = {^eM.'^,\^\< ^J. ■ 

Now, we bound the conditional variance (condition iii. in Theorem [3]) . 

Lemma 4 Assume [W\) and fix to,M > 0. Then for every s,t £ [to^T] and x,y G 
[-M,M]'^ 

Var{u{t, x) \u{s, y)) > CA{{t, x); (s, y)) 



where A is the metric given by i31\) . 

Proof: We wih use the fohowing formula: if (U, V) is a centered Gaussian vector, then 



(Puy - i^^u - (^v?_){{(^u + (yy? - ply) 

4C7? 



Var{U, V) = Y^ (32) 



'V 
where pfjy = E(C/ - Vf, a^ = EC/^, a^ = EV^. Denote by 

Plx,s,y = '^Wit,x) -u{s,y)f , a^{t,x) = Bu{t,xf, aly = Eu{s,yf. 
It suffices to show that 

iplx,s,y - {^t,x - Crs,yf){{(Tt,x + (ys,yf - Pt,x,s,y) > cA((t,x); {s,y)) 

for every s,t S [ioi 2^] and x, y G [— M, M] . By Theorem [2] the second factor in the left-hand 
side above is bounded below by a constant. So it remains to check that 

iplx,s,y - K,x - (ys,yf) > cA((t,x); (s,y)) 
but this has been done in pTC], proof of Proposition 3.2. (see also [U], proof of Lemma 4.3). 
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Remark 7 Using the previous result we can give a bound on the joint density pt,x,s,y of 
the vector {u{t,x),u{s,y)). Actually, one can show that for every t € [to,T] and x,y G 
[— M, M]'^ we have the inequality 



Pt,s,x,y{zi, Z2) < CiA((t,x);(s,y)) 2 exp f - 



\{it,xy,is,y)) 



for every zi,Z2 G [—N,N]^, where A is the metric defined by Ii31\) . It suffices to follow the 
lines of Proposition 3. 2 in |J5| /. 

We can state now the main result of this section. 

Theorem 4 Let us consider I, J non-trivial compacts sets in \to,T] and [—M,M] respec- 
tively. Fix N > and let u be the solution to the system ([I]j. Then for every Borel set A 
contained in [—N,N] it holds that 

C-'^ Capk-^ < P {u{I X J) n A / 0) < CHk-^iA) 

with 

2{d + l) 

J = k 



2H + 1-/3' 

Proof: The proof is a consequence of Theorem [3] and of the preceding two lemmas. 

■ 

Remark 8 

• Of course, for H close to 2, our result recovers the findings in J^. 

• it is also possible to give some results concerning the probability that, for fixed t, x, the 

sets u{{t} X J) and u{I x {x}) (as before I, J non-trivial compacts sets in [to,T] and 
in [—M,M]'^ respectively) to hit a given Borel set A contained in [—N,N]''. Actually, 
by routine arguments we will have 

C'^ Cap, 2d (A) < P (u({t} xJ)r\A^m< Cn. 2d (A) 

" 2H+l-fJ ''-2H+1-I3 

and 

C-^Cap, 2 (AXP (u(I X {x\) n A 7^ 01< Cn. 2 (A). 
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